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Abstract 

The cell membrane is inherently asymmetric and heterogeneous in its composition, a fea- 
ture that is crucial for its function. Using atomistic molecular dynamics simulations, the phys- 
ical properties of a 3-component asymmetric mixed lipid bilayer system comprising of an un- 
saturated POPC (palmitoyl-oleoyl-phosphatidyl-choline), a saturated SM (sphingomyelin) and 
cholesterol are investigated. In these simulations, the initial stages of liquid ordered, l Q , do- 
main formation are observed and such domains are found to be highly enriched in cholesterol 
and SM. The current simulations also suggest that the cholesterol molecules may partition into 
these SM-dominated regions in the ratio of 3 : 1 when compared to POPC -dominated regions. 
SM molecules exhibit a measurable tilt and long range tilt correlations are observed within the 
l a domain as a consequence of the asymmetry of the bilayer, with implications to local mem- 
brane deformation and budding. Tagged particle diffusion for SM and cholesterol molecules, 
which reflects spatial variations in the physical environment encountered by the tagged parti- 
cle, is computed and compared with recent experimental results obtained from high resolution 
microscopy. 
*To whom correspondence should be addressed 
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Introduction 



The cell membrane is characterized by both lateral and transverse lipid heterogeneity, an aspect of 
significant functional consequence. 1 Transverse lipid heterogeneity is maintained actively by the 
cell, making the cell bilayer intrinsically asymmetric. Lateral lipid heterogeneities called 'rafts',^ 
which are ternary mixtures of sphingomyelin (SM), phosphotidyl Choline (PC) and cholesterol 
(Choi) molecules have been implicated in a variety of cellular processes including signaling and en- 
docytosis, even though the nature of these functional cellular rafts is still a matter of contention.®^ 
Asymmetry in bilayers can arise both in terms of difference in constituent lipids or in number 
of lipid molecules in both leaflets. In spite of this obvious lateral/transverse compositional hetero- 
geneity, except for a few seminal studies,^ most in- vitro investigations of multicomponent artificial 
membranes'^^have been done on symmetric bilayers. Further, most of the atomistic simulations 
of model cell membrane mimics, have been carried out on systems which either have lateral het- 
erogeneity^or transverse asymmetryl-^^but rarely both. There have however been a few studies 
using coarse-grained simulations J^lHI and continuum Landau theories^^that address the role of 
inter-bilayer coupling in equilibrium physical properties and domain growth. 

While early simulations of model membranes consisted only of a single-component PC bilay- 
ers Ej}|23] i a (- er simulations have incorporated more than one lipid component, in particular choles- 
terol. 24 Following the 'raft proposal' of the importance of SM molecules in raft formation, an 
increasing number of simulations with SM molecules'^ have been carried out. These simula- 
tions include detailed comparisons between mixtures of SM and Choi and PC and CholM 2 ^^ 
Finally, simulations on asymmetric bilayers have also been reported. A simulation study by Bhide 
et al.^ was performed on systems consisting of SM and Choi molecules in upper leaflet and SOPS 
and Choi molecules in the lower leaflet, with equal number of phospholipids in both the leaflets. 
Comparisons between each leaflet of the asymmetric bilayer with corresponding simulations of 
symmetric bilayers showed no significant differences in their physical properties. 

In this paper, the equilibrium properties of a multicomponent asymmetric (both lateral and 
transverse) bilayer using atomistic molecular dynamics (MD) simulations are studied. Specifically, 
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the physical properties of a 3-component asymmetric lipid bilayer comprising of an unsaturated 
POPC, a saturated SM and Choi molecules, which exhibits lateral compositional heterogeneities 
in the form of liquid ordered (l ) - liquid disordered (l c i) domains, thought to represent the charac- 
teristic lipid composition of rafts on the cell membrane. In addition to studying a variety of order 
parameters, their correlations and spatial distribution, transport properties of the component lipid 
molecules are also studied. This study suggests that the presence of lateral heterogeneities in the 
bilayer can potentially affect the molecular diffusion at short time scales, which may be of rele- 
vance to molecular diffusion at the cell surface. However lipid based lateral heterogeneities on the 
cell surface are likely to be small and below optical resolution; it is only recently that advances in 
high resolution single -particle tracking (SPT)^^2 allow one to measure subtle changes in diffusion 
characteristics as molecules move across the heterogeneous cell surface. In addition, the transverse 
asymmetry would suggest that the diffusion characteristics of the molecule is different in the two 
leaflets of the bilayer. Correlations between the local physical environment and the tagged particle 
diffusion are also investigated in this paper. 

Materials and Methods 

Two bilayer model systems have been simulated : (A) a symmetric three-component bilayer con- 
sisting of 170 POPC, 171 palmitoyl-sphingomyelin (SM, saturated lipid) and 171 Choi (ratio of 
1 : 1 : 1) in each leaflet and (B) an asymmetric bilayer with a composition of 170 POPC, 168 SM 
and 171 Choi (ratio of roughly 1 : 1 : 1) in the upper leaflet and 256 POPC and 256 Choi (ratio 
1 : l)in the lower leaflet. Two control systems have also been simulated, which are (C) a symmetric 
one-component bilayer made up of 128 pure POPC and (D) a symmetric two-component bilayer 
made up of 64 POPC and 64 Choi (ratio 1 : 1 in both leaflets). 

The force field parameters for POPC were taken from the previously validated united-atom 
description (Tielman and Berendsen 3 ^), whereas for SM and Choi, parameters from works of 
Niemela et al.^ were used. Water was simulated by the simple point charge (SPC) model^ and 
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PACKMOL^was used to generate the initial configurations of all the bilayer systems and hydrated 
with water in the ratio of 1 (lipid) : 32 (water). 

All MD simulations were performed at a temperature T = 296K, which lies between the main 
transition temperatures of POPC (T m = -2.9 ± l.3°Cf® and SM (T m w 41.4°C) 38 and right in 
the middle of the l -ld phase coexistence region of the symmetric ternary system. 39 GROMACS 
(http://www.gromacs.org) 40 software was used to integrate the equations of motion with a time 
step of 2fs. To avoid bad contacts arising from steric constraints during initialization, all bilayer 
systems were subjected to steepest descent minimization initially. The systems were then simulated 
for 50 ps in the NVT ensemble using a Langevin thermostat. Subsequently, each system was 
simulated in the NPT ensemble (T = 296K, P = latm) using a Berendsen thermostat and semi- 
isotropic pressure coupling with compressibility 4.5 x 10~ 5 bar -1 for 220 ns (for the asymmetric 
ternary system B) and for 100 ns (for symmetric ternary system A) and for 100 ns (for the single and 
two -component systems C and D). Lipid system parameters such as deuterium order parameter 
and mean thickness were monitored throughout the simulations to ensure that the system is well 
equilibrated. All the reported results are computed over last 20ns unless otherwise stated. 

Note that the two model systems A and B have been simulated starting from two sets of initial 
conditions : (i) where the components in each leaflet are mixed or (ii) where the ternary compo- 
nents are completely phase segregated. 



Results 

To make sure that a stable, surface tension-less asymmetric bilayer is being simulated, the forces, 
torques and surface tension of the bilayer are computed from the local stress tensor Oij(x,y,z) = 
v £a f^ffi where ff* is the i th component of the force on the a th -particle due to all other particles 
within a coarse-grained volume v = [0. 1 nm] 3 . The net force F, = J dkO^dv and its first moment 
(related to the torque) = J (d/Onx^ — dia^Xi) dv for the bilayer were computed and both force 
and torque balance - F\ = 0.29 ± 3.34 , F 2 = -2.56 ± 2.09, F 3 = -2. 16 ± 3.2, in units of nN and 
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Mi 2 = 0.77 ± 0.428, M 13 = -8.83 ±2.04, M 23 = -13.303 ± 1.67, in units of nN-nm were achieved 
suggesting a mechanically stable asymmetric bilayer. These are comparable to the corresponding 
values for the symmetric bilayers (see Supplementary Information). The membrane surface tension 
is calculated as y = J n(z)dz, integrated over the width of the bilayer, where n(z) is the lateral 
pressure, given by n(z) = \ {p xx {z) + o yy (z)) — a zz (z) resulting in a value of 7= —0.0018 ± 
0.0301 bar-nm, essentially a 'zero' surface tension bilayer. More details of time dependence of net 
forces, their moments and surface tension are provided in Supplementary Figure S2. A snapshot 
of the equilibrium bilayer configuration of the asymmetric membrane is shown in Figure 1(a). In 
Figure 1(b), the pressure profiles of symmetric and asymmetric bilayers are shown for comparison. 
In contrast to the asymmetric bilayer, the symmetric bilayer exhibits a symmetric pressure profile 
about the midplane. On the other hand the pressure profile of the asymmetric bilayer shows larger 
spatial oscillations. 

Lipid composition 

At the temperature and overall lipid composition under consideration, the homogeneous mixed 
phase of both the asymmetric and symmetric bilayers is unstable and exhibits definite features of 
phase separation between POPC-rich and SM-rich domains (Supplementary Figures, S3, S4); in 
the asymmetric bilayer phase segregation occurs in the upper leaflet alone, while the composition 
in the lower leaflet remains homogeneous. To show that the bilayer is undergoing phase separation 
towards a final equilibrium phase segregated configuration, the theory of dynamical coarsening,^ 
which deals with the study of the dynamics of domain formation starting from a complete disor- 
dered phase is used. In simulations starting from an initial mixed state, it was found that the system 
quickly phase segregates; at early times the domain sizes are small and grow over the simulation 
time, 220 ns. It has been shown earlier that complete phase segregation occurs only over a very 
long time, at least on a microsecond time scaleJ^EI To demonstrate dynamics of domain forma- 
tion or coarsening, we have computed (i) the probability distribution P((j)) of the order parameter 
= p™+pp°^ • at different times (where p SMjmPC is the local density of SM/POPC), and (ii) the 
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energy density e = E /V, where the energy E /d 2 r(V0) 2 , which measures the time dependence 
of the amount of interface separating the two phases. The results are shown in Figure 2. Initially 
P((j>) is peaked at = (mixed state) and subsequently evolves into a distribution with two peaks 
at ±1 which gets progressively sharper with time. Dynamical scaling in the coarsening regime^ 
implies that e ~ t~ l ^ z , where z = 3, since the order parameter is conserved. Our simulations show 
that l/z = 0.30 ± 0.15, consistent with this scaling prediction (Figure 2 (a)). 

The domains are observed to be larger in the symmetric bilayers compared to the asymmetric 
bilayer, when measured over the total simulation time scale (220 ns). The interfacial fluctuations 
of the / -domains in the asymmetric bilayer are larger than the symmetric bilayer, suggesting, in 
line with earlier simulations, 45 that the interfacial tension between the domains is higher in the 
symmetric bilayer. Both the larger interfacial tension and the enhanced correlations between the 
domains in the two leaflets of the symmetric bilayer (seen in Supplementary Figure, S4), are a 
result of a transbilayer coupling. This transbilayer correlation between domains can be understood 
within a Landau description of both equilibrium and dynamics, in which the relative concentration 
of lipids in one leaflet acts as a "field" for the concentration in the other leaflet. 20 Local equilib- 
rium relates the domain growth velocity to the interfacial tension, 43 this suggests that the domains 
coarsen faster in the symmetric bilayer, consistent with our simulations. 

To estimate the relative partitioning of cholesterol in the POPC and the SM rich regions in 
the upper leaflet, the joint probability distribution of finding a given concentration of SM with 
Choi in an xy region (similarly, POPC with Choi) is calculated and shown in Figure 2 (b) and 2 (c). 
These joint probabilities show that the cholesterol concentration completely correlates with the SM 
concentration (i.e., Choi is low (high) when the SM concentration is low (high)) and completely 
anticorrelates with the POPC concentration (i.e., Choi is low (high) when POPC is high (low)). 
The joint probability densities plotted in Figures 2(b) and (c) strongly suggest that cholesterol 
preferentially partitions in the SM-rich phase three times more than in POPC -rich region (more 
precisely 2.97 : 1) 

The segregation of chemical composition is accompanied by changes in physical characteristics 
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of the segregated molecules. The saturated lipid tails of SM that are in the SM-enriched domains 
are found to be more rigid than both the SM and POPC molecules in the POPC-rich domains. The 
deuterium order parameter (S) describing the rigidity of lipid tails is computed for SM molecules 
in both the leaflets and is shown in Figure 4. The probability distribution, P(S), of SM gives 
indication of the location of SM in either SM-rich (higher value) or POPC-rich domains (lower 
value). As seen in Figure 3(a), P(S) has a bimodal distribution. The distinction of the deuterium 
order parameter in the two regions indicate that the SM- and POPC-rich regions may be identified 
with Z -like and Zj-like phases, respectively. From Figure 3(a), it can be seen that the probability 
distribution of deuterium order parameters in the asymmetric bilayers are comparable to that of 
symmetric bilayer. The spatial variation of the deuterium order parameter in the x—y plane for both 
asymmetric and symmetric bilayers is shown in Supplementary Figure S7 and it can be noted that 
the segregation of chemical composition is naturally accompanied by an / -/,./-like phase separation 
in the bilayer membrane. From this spatial distribution, it can be seen that the size of / c -like 
domain size (enriched in SM) is significantly larger in the symmetric bilayer compared to that of 
asymmetric bilayer, over the time scale measured. 

The probability distribution and the spatial profile of the local bilayer thickness is also com- 
puted and shown in Figure 3(b) and Supplementary Figure S7 respectively. The SM-enriched 
/ -like domains have a larger bilayer thickness compared to the POPC- enriched Z^-like domains. 
As with deuterium order parameter, the distribution of bilayer thickness for both the symmetric and 
asymmetric bilayers is bimodal, corresponding to the /^-like and Z -like domains. It is significant 
that the difference in the bilayer thickness between the / -like and Z^-like phases in the symmetric 
bilayer is consistent with recent AFM studies, 46 and is more than twice compared to asymmetric 
bilayer. The spatial profile of the thickness shown in Supplementary Figure S7 is an accurate mea- 
sure of the domain size and consistent with the deuterium order parameter results, the domain size 
in the symmetric bilayer is larger than that of asymmetric bilayer over the time scale measured. 

By computing the joint probability of finding SM molecules in the upper and lower leaflets of 
the symmetric bilayer, it can be seen (Supplementary Figure S5) that there is a clear registry of 
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SM-rich (l -\ike) domains across the leaflets, which may be responsible for observed increase in 
SM-rich domain sizes in symmetric bilayers. 

Lipid splay and tilt 

To quantify the relative packing of lipid chains, the amount of splay between the two lipid tails 
is calculated, in addition to deuterium order parameter computed above. A tail vector is defined 
as a vector originating from the carbon of the carbonyl group and pointing to the terminal methyl 
carbon of a lipid tail. The splay angle is the angle between two such tail vectors of each lipid. 
Simulation results suggest that the extent of splay in SM lipid tails is significantly smaller in l - 
like domain (17°) compared to POPC-rich domain (40°) as seen in Supplementary Figure S6. 
In addition, the extent of lipid splay in the symmetric and asymmetric bilayers is comparable. 
The spatial heterogeneity of the splay in the asymmetric bilayer reflects the differences in lipid 
composition (Supplementary Figure, S8). Taken together with the deuterium order parameter data, 
the lipid splay results strongly suggest that the packing fraction of lipids in the SM-rich phase is 
higher than that in the POPC-rich phase. 

Tilt angle of the lipid tail chain is defined as the orientation of the mean tail vector of the lipid 
with respect to the local outward normal to the membrane, and is described by two angles (0,0), 
the polar and azimuthal angles, respectively. The angle measures the orientation of the 2d tilt 
vector, the projection of the tail vector onto the tangent plane, with the jc-axis. An accurate deter- 
mination of the tilt angles of the component lipids is quite involved, since, over short length scales, 
the local normal to the membrane fluctuates due to molecular protrusion effects. To compute lo- 
cal average tilt, a coarse-graining scale is chosen, which should be more than the protrusion scale 
and less than the tilt correlation length (which in turn should of course be smaller than the size of 
the /„-like domain). Over this coarse-grained scale, a membrane normal is considered to be along 
the z-axis. A convenient choice of coarse graining scale is around a 1 ran 2 (which encompasses 
3 lipids on an average), for which statistically reliable results can be obtained. For instance, for 
POPC-only bilayer, the probability distribution of the coarse-grained angle, P(6), is peaked about 



8 



zero, while the distribution P((j)) is uniform, consistent with the known fact that POPC does not 
exhibit a tilt at this temperature. 

Following such a procedure for computing local tilt angles for the ternary membrane is prob- 
lematic, since over the timescale of the MD simulation, the size of the Z -like domains are small, 
making it difficult to obtain an unambiguous determination of tilt and its correlations. To address 
this issue, simulation of a ternary bilayer starting from a completely phase segregated configura- 
tion with cholesterol distributed, between SM-rich and POPC-rich domains, in the ratio of 3:1 (as 
described in Materials and Methods) is carried out. In this case the Z -like domain (SM-rich) size 
is around half the system size and can provide good statistics. 

The tilt angle distribution of POPC lipids in the symmetric and asymmetric ternary bilayers 
has a similar distribution as in the case of POPC-only bilayer as seen in Supplementary Figure 
S9. On the other hand, the tilt angle distribution for SM shows an interesting trend. While SM 
in the symmetric bilayer shows no evidence of tilt (P(6) is peaked at around 0° and P((j)) is 
uniform), SM in the SM-rich (Z -like) domain of the asymmetric bilayer has a nonzero tilt of around 
2.016 ±0.69°. This small tilt of SM in the asymmetric bilayer is consistent with the decrease in 
bilayer thickness of the asymmetric membrane compared to that of symmetric bilayer (see Figure 
3(b)). Tilt angle correlations defined as C(r) = (0(r)0(O)) are computed for POPC and SM lipid 
molecules and the results are shown in Figure 4. The correlation functions of the tilt of POPC and 
SM (symmetric bilayer) decay exponentially to zero, consistent with the above findings. However, 
C(r) for SM in the asymmetric bilayer decays exponentially to a non zero value, signaling long 
range order^in SM molecules in the asymmetric bilayer. This suggests that the asymmetric nature 
of the bilayer can potentially generate a tilt ordering of SM molecules in the SM-rich (/ -like) 
domain. 

Previous simulations^ 1 2 I 45 l 48 l and experiments^lH^lhave shown that the thickness of l Q domains 
is larger than that of lj domain. This difference in thickness between l and l c \ domains gives rise 
to line tension along the domain boundary.^ The lipid tails of the l D domain can be exposed to 
solvent as a result of such mismatch, which is energetically unfavourable. One of the ways to 
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mitigate such mismatch is for the thicker l D domain to undergo a small tilt such that the head 
groups of the two domains can be at the same height. The observed small tilt in the simulations of 
asymmetric bilayers supports this hypothesis. Within a Landau theory, the asymmetric bilayer can 
be thought of as being subjected to a transverse compression, which would naturally lead to a tilt 
when the lipids are stretched out (as they are in the Zo domain) ElES 

The existence of a finite tilt and its correlation over long scales, if verified experimentally,^ 
could have important consequences for membrane deformation and budding The tilt vector nat- 
urally couples to the local curvature tensor of the membrane, giving rise to anisotropic bending 
stresses. If the constituent molecules are chiral (as they usually are), then there are additional 
bending stresses coming from chiral couplings of the tilt and curvature. If strong enough, these 
bending stresses can induce membrane deformation giving rise to spherical buds or cylindrical 
tubules.EEH 

Tagged particle diffusion 

Even though the membrane composition forms stable domains at equilibrium, the individual com- 
ponent molecules can traverse across domains. The equilibrium dynamics of certain 'tagged' com- 
ponent lipids can be monitored by measuring their diffusion coefficients and correlating them with 
the physical and chemical heterogeneity across the membrane. The mean square displacement 
(MSD) is defined as (<5r,-(?) 2 ) of a tagged particle, where 5r,-(f) = r,-(f) — r,-(0) is the displacement 
of tagged i th lipid of a given species at time t from its position at t = 0. While computing the MSD, 
the location of the tagged lipid, whether it is in the POPC-rich (Z^-like) or SM-rich (Z -like) do- 
main, is monitored by computing the instantaneous deuterium order parameter S of the tagged lipid 
and the local bilayer thickness. The diffusion analysis is show in Figure 5. It can be seen that the 
diffusion coefficient of both SM and Choi molecules, given by the saturation value of (8rf(t))/4t, 
is different depending on whether the tagged lipid molecule is in the Z -like and Z^-like domain as 
has been reported in other simulation studies.^ 

For tagged SM or Choi molecules that initially lie in the POPC-rich (Z c ;-like) domain, the local 
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diffusion coefficient crosses over from an early time high value (Do) to a late time low value (£>«,). 
For each of the tagged components, a first-passage time, defined as the first time that a tagged 
molecule residing in a domain moves out of it, is computed. The crossover time % is then obtained 
from such a computed first-passage time. The MSD is found to obey a crossover scaling relation, 

(8r 2 )=4D tF(t/T) (1) 

where the nonlinear scaling function F has the asymptotic form, 

F(t/x) = 1 for f/r<l 

= — for f/r>l (2) 
Do 1 

The data collapse shown in Figures 5(a) and 5(b) demonstrate this crossover scaling for SM and 
Choi molecules, respectively. Experimentally the crossover time scale % can be obtained from the 
value of t at which the instantaneous D(t) = (Dq+Doo) /2 , rather than the first-passage time.The 
values of the diffusion coefficients for SM molecules (Fig. 5(a)) are in agreement with the ex- 
perimental findings, ^S159| w hii e fa e values for Choi molecules (Fig. 5(b)) are slightly lower than 
those reported in other simulations.^ This discrepancy can be attributed to the differences in lipid 
composition of the ternary bilayer used in both the simulations. The difference in the tagged par- 
ticle diffusion coefficient between the two domains can be attributed to changes in local viscosity 
(t]), moment-of-inertia (/) of the particle and to changes in the local density correlations (given 
by the local partial structure factors, S a p (q)) with neighbouring molecules that the tagged particle 
experiences as it traverses across the domain. The change in the deuterium order parameter S is 
used to track the changes in moment-of-inertia of the tagged molecules. For SM molecules, this 
shows a crossover similar to the crossover diffusion, with the value of S being low (high) when 
the diffusion coefficient is high (low) as seen in Figure 5(a). On the other hand cholesterol being a 
rigid molecule is unlikely to undergo any conformational change, resulting change in moment-of- 
inertia, as it traverses across the domains. Thus a substantial change in the diffusion coefficient of 
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Choi molecules can be attributed to the changes in viscosity and local density correlations, S a p (q), 
arising from changes in the local environment. 

In the context of SPT experiments, if the time scale over which the tagged particle is tracked 
is large enough so that the particle crosses and recrosses the domains, then one would measure 
the MSD of 8ri(t) = J dt' {n{t' +t) -r* (*')). For the control pure POPC symmetric bilayer, it 
is seen that the trajectories are purely diffusive ((5r;(f) ) t) in nature. For the ternary bilayer, 
however, this MSD would exhibit deviations from true diffusion, which can be characterized by 
(8ri(t) 2 ) oc t a . Figure 5c shows the MSD for the tagged lipids in both the symmetric and asym- 
metric bilayers and the corresponding values of the exponent a. The value of a obtained for SM in 
the asymmetric bilayer (a ~ 0.39) is consistent with the experimental value obtained by analyzing 
recent SPT of labeled SM on the plasma membrane of epithelial cells, a«0.3. 33 The interbilayer 
coupling in the symmetric ternary system, makes the tagged SM movement slower, as seen by the 
lower value of a ~ 0.35 (see Figure 5(c)). The tagged particle dynamics of the other lipids in 
the asymmetric bilayer, viz., POPC and cholesterol, show interesting differences between the two 
leaflets - typically, the upper leaflet lipids show a smaller a (a = 0.45,0.41, respectively), than 
the lower leaflet lipids (a = 0.85,0.89, respectively). 

Conclusion 

The cell membrane exhibits both lateral and transverse heterogeneity. In this paper, the equilibrium 
properties of a ternary component asymmetric bilayer membrane system at l -lj phase coexistence 
are studied using an atomistic MD simulation over a time scale of 220 ns. The asymmetric bilayer 
considered in this study is composed of POPC, SM and cholesterol in the ratio of 1 : 1 : 1 in the upper 
leaflet and POPC and cholesterol in the lower leaflet. The two significant results from this study 
are: (i) cholesterol prefers to be associated with SM-rich domains (Z -like) three times more that 
POPC-rich domains (/^-like) and (ii) the saturated lipid SM in the / -like domain exhibits long- 
range tilt correlations purely as a result of the asymmetry in the bilayer composition (in contrast, 
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the SM lipid molecules in the symmetric bilayer show no such tilt). 

This bilayer- asymmetry induced lipid tilt in the Z -like domain has important implications to 
local membrane deformation and hence membrane budding and endocytosis. The existence of 
a lipid tilt expressed over large scales provides a natural coupling to the local curvature tensor, 
and results in anisotropic bending stresses at the membrane.^ Moreover, if the constituent lipids 
are chiral (as they are in 'raft' -lipids), then there would be additional bending stresses serving to 
deform the membrane locally. When the strength of these couplings are large enough, they can 
induce the local formation of spherical buds or cylindrical tubules. Recent FRET-based studies 
of the organization of lipid tethered proteins on the outer surface of living cells, such as GPT 
anchored proteins, have shown that they form cholesterol sensitive nanoclusters mediated by the 
activity of cortical actin.^1 These studies imply that there must exist a molecular linkage between 
the outer leaflet GPTanchored proteins and cortical actin. The current study will form the basis 
for further investigations on possible transbilayer interactions between GPTanchored proteins, SM 
and cholesterol, with specific saturated, long chain lipids at the inner leaflet that have potential 
interactions with actin or actin remodeling proteins. 
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(a) (b) 



Figure 1: (a) Snapshot (side view) at the end of the simulation of the ternary system POPC (gray), 
SM (orange), Choi (yellow) forming the stable asymmetric bilayer in water (cyan), (b) Lateral 
pressure profiles n(z) as a function of position across the bilayer, for the ternary asymmetric bilayer 
(blue), and the ternary symmetric bilayer (red). 
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Figure 2: (a) Probability distribution P((j)) (normalized) of the order parameter = 
where p SM and p POPC are the concentrations of SM and POPC, respectively. Its time dependence from 
t = — 200 ns, shows that the system, initially prepared in the mixed state (0 = 0), has coarsened 
into SM rich (0 = 1) and POPC rich (0 = — 1) domains separated by sharp interfaces. The last 
panel shows the dependence of the energy density £ with time. The data (red filled circles) suggests 
that £ goes as a power-law, £ ~ f--3±0.i5 (fi^ed line), with an exponent consistent with dynamical 
scaling, z = 3 (see text), (b) Joint probability distribution (color bar) of the concentration (in units 
of number/nm 2 ) of (b) SM and Choi and (c) POPC and Choi, coarse-grained over [1.73 nm] 2 in 
the upper leaflet and averaged over 20 ns. Red shows the highest joint probability and blue the 
lowest. Note the strong correlation (anticorrelation) between SM-Chol (POPC-Chol), respectively. 
The figures clearly show an enrichment of cholesterol in SM-enriched domains by a ratio 3 : 1 
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Figure 3: (a) Probability distribution of the deuterium order parameter S in the asymmetric (blue) 
and symmetric (red) bilayers. This is then coarse-grained over a spatial scale of 1.56 nm and time 
scale of 20 ns. Here S of the selected carbons (C5-C7) of POPC and SM is displayed, after suitable 
binning in the xy-plane. with the lower branch 0.1 —0.26 corresponding to the low partitioning 
of SM in the phase and the higher branch 0.35 — 0.45 corresponding to its enrichment in the l 
phase, (b) Probability distribution of the bilayer thickness d in the asymmetric (blue) and symmet- 
ric (red) bilayers. The two distinct peaks in the asymmetric bilayer at 4nm and 4.3 nm indicates 
the coexistence of the 1^ and l phases. The local bilayer thickness is constructed by via x-y grids 
for the upper and lower leaflet of the bilayer, with grid size = 1.95 nm. An average over the z- 
coordinates of the phosphorus atom in the head groups of POPC and SM within each x-y grid is 
then performed. The difference between the z-coordinates corresponding to the same x-y grid of 
the two leaflets gives the local bilayer thickness. 
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Figure 4: Probability distribution of the coarse-grained tilt angles 6 and of SM (see text) in (a, 
b) the asymmetric bilayer, and (c, d) symmetric bilayer, respectively. These data are collected from 
equilibrium configurations that exhibit complete phase segregation, where the domain size is about 
half the system size. In the asymmetric ternary bilayer, SM exhibits a measurable tilt, as evidenced 
from the distributions of 6 and 0. In comparison, SM in the symmetric bilayer shows an absence 
of tilt - the distribution of is peaked at a value close to 0, while the distribution of is uniform, 
(e) The tilt correlation function C(r) = (0(r)0(O)}, normalized to its value at r = 0, shows an 
exponential decay to zero for POPC (asymmetric bilayer) and SM (symmetric bilayer). However 
SM in the asymmetric bilayer shows an exponential decay to a nonzero value, demonstrating long 
range tilt correlations. The correlation length, defined by the scale at which C(r) decreases to \ je 
of its value at r = 0, can be easily read out, t; SM (sym) = 0.35 nm, t; SM (asy) = 0.5 nm and t, m?c {asy) = 
0.3 nm. 
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Figure 5: Mean square displacement of tagged particles as a function of time collected at high 
resolution over short time scales. Tagged particle diffusion of (a) SM and (b) Choi in the asym- 
metric bilayer- statistics collected over 6 tagged particles, each starting from the lj domain, shows 
a data collapse onto a nonlinear crossover scaling curve. The crossover from high to low dif- 
fusion occurs at a time z, computed as the first-passage time of the tagged particle. For SM, 
we find that D = 3.49/jm 2 s _1 and = 2.69/im 2 s -1 , while for Choi, D = 4.23 J um 2 s~ 1 and 
Doo = 1.36jUm 2 s _1 . Note that in (a) as SM moves across domains, the deuterium order parameter 
S of tagged SM shows a similar crossover, going from ~ 0.22 (Z^) to ~ 0.42 (/ ). (c) MSD ver- 
sus time with fits to (5r ; (f) 2 ) oc t a , for SM in the asymmetric (a = 0.39) and symmetric bilayers 
(a = 0.35), and POPC and Choi in the upper leaflet of the asymmetric bilayer (a = 0.45 and 
0.41, respectively). Note that the a exponent for SM in the symmetric bilayer is smaller than the 
asymmetric bilayer, due to the transbilayer coupling. Likewise, the MSD for POPC and Choi in 
the symmetric ternary bilayer are shown in Supplementary Figure S10, for comparison. 
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Supplementary Information 
A. Polley et al. 

Tabulation of force, torque and tension in model bilayers at equilibrium 

1. Components of Force 

POPC (symmetric) 

F x = 0.83 ±1.66nN,F 2 = 0.27 ±1.14nN,F 3 = -0.16±1.25nN 
POPC-CHOL (symmetric) 

F x = - 1 .01 ± 1 .36 nN, F 2 = -0.32 ± 2.72 nN, F 3 = 0. 18 ± 0.62 nN 

2. Components of the moment of force 

POPC (symmetric) 

M12 = 8.26 ± 0. 102 nN-nm, M i3 = 10.74 ± 0.066 nN-nm, M 23 = 9.33 ± 0.094 nN-nm 
POPC-CHOL (symmetric) 

M12 = 5.94 ± 0.082 nN-nm, M i3 = 8.882 ± 0.23 nN-nm, M 23 = 10.757 ± 0.42 nN-nm 

3. Surface Tension 

POPC (symmetric) 

7 = 0.0204 ± 0.0357 bar-nm 

POPC-CHOL (symmetric) 
7 = 0.0054 ± 0.0099 bar-nm 

The values of the components of force, moment of force and surface tension in these sym- 
metric 1 and 2 component systems are comparable to those of the asymmetric ternary system 
given in the main text. This implies that the asymmetric ternary bilayer is as mechanically 
stable as the symmetric systems. 
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Supplementary Figures 
A. Polley et al. 




Figure 1: Lateral pressure profile of the symmetric bilayer made of (a) pure POPC (blue) and (b) 
POPC + Choi, with ratio 1 : 1 (red). 
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Time (ns) Time (ns) Time (ns) 

Figure 2: (a) Time dependence of the components of the force (defined in main text), where each 
data point is an average over 10 ns . (b) Time dependence of the components of the torque (defined 
in main text) where each data point is averaged over 10 ns. The observation that at late times these 
components fluctuate about zero, and that these fluctuations are not large, implies that in the mean 
there is force and torque balance, (c) The time dependence of the surface tension (evaluated from 
the components of the stress tensor) in the upper (y«) and lower (7/) leaflets of the membrane and 
the total surface tension (7). The observation that both y u and 7/ are small and fluctuate about 
zero, implies that in the mean there is no imbalance of surface tension in the two leaflets. The 
observation that the total 7 is small and fluctuates about zero implies that in the mean the bilayer 
membrane is tensionless. 
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Figure 3: Spatial variation of lipid concentration in the asymmetric ternary bilayer system : (a) 
POPC, (b) PSM, (c) Choi in the upper leaflet, and (d) POPC and (e) Choi in the lower leaflet. The 
LUT bars denote the fraction of the lipid species within an area [1. 56 nm] 2 . 
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Figure 4: Spatial variation of lipid concentration in the symmetric ternary bilayer system : (a) 
POPC, (b) PSM, (c) Choi in the upper leaflet, and (d) POPC, (e) PSM and (f) Choi in the lower 
leaflet. The LUT bars denote the fraction of the lipid species within an area [1. 56 nm} 2 . Note from 
(b) that the domain sizes are larger than in the asymmetric bilayer. Further, from (b) and (e), we 
note that there is significant registry of the domains across the bilayer (see Fig. 5). 
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Figure 5: Joint probability distribution (color bar) of the concentration (in units of number / nm 1 ) 
of SM of upper leaflet and SM at lower leaflet of symmetric ternary bilayer coarse-grained over 
[1.73 nm] and averaged over 20 ns. Red shows the highest joint probability and blue the lowest. 
Note the strong correlation between SM-SM in the upper and lower leaflets implies significant 
bilayer registry in the SM-enriched domain across the two leaflets. 
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Figure 6: [A] Probability distribution of the deuterium order parameter S of POPC in the symmetric 
bilayer of (a) POPC and (b) POPC + CHOL (with ratio 1:1). [B] Probability distribution of 
splay of POPC in the symmetric bilayer of (a) POPC and (b) POPC-CHOL (with ratio 1:1). [C] 
Probability distribution of splay of (a) POPC and (b) PSM in the upper leaflet and (c) POPC in the 
lower leaflet of the asymmetric ternary bilayer. 
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Figure 7: Spatial variation of deuterium order parameter S(x,y) and bilayer thickness d(x,y) of 
(a-b) asymmetric bilayer and (c-d) symmetric bilayer of the ternary system. In both cases, the 
domains coarsen over time; note however that the domain are larger in the symmetric bilayer 
compared to the asymmetric bilayer, evaluated over the same time. In addition, note that in the 
case of the symmetric bilayer, the domain registry across the bilayer is significant, consistent with 
Fig. 5. 
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Figure 9: Probability distribution of the tilt angles (a) 6 and (b) for POPC in the asymmetric 
ternary bilayer. Note that the distribution of is peaked about 0, while the distribution of is 
uniform, consistent with the known fact that POPC does not exhibit a tilt at these temperatures. 
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Figure 10: MSD versus time for tagged POPC and Choi in the symmetric ternary bilayer. Fits 
to (8rj(t) 2 ) oc t a gives a = 0.35 (POPC) and a = 0.32 (Choi). Note that the values of a are 
consistently lower than those in the asymmetric ternary bilayer. 
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